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Abstract 

The best available heat-transfer and friction corre- 
lations for single-phase turbulent flow of hydrogen 
through tubes with large ratios of wall to bulk tempera- 
tures, are applied to the cooling passages Of nuclear roc- 
ket nozzles. The calculated pressure drop and temper- 
ature rise of the coolant are compared to the measured 
values for nuclear test runs of Pewee, NERVA, and 
Phoebus 2A nozzles. Calculated local gas-side wall tem- 
peratures are compared with wall temperatures indi- 
cated by melting of selected braze alloys from NERVA 
nuclear tests. The use of several different heat -transfer 
correlations on the hot-gas side and their effect on the 
predicted wall temperature is also discussed. 

I. Introduction 

An effective method of predicting heat-transfer and 
friction coefficients in the coolant passages of regener- 
atively cooled rocket nozzles is essential to the opti- 
mization of any nozzle design. In nuclear rocket noz- 
zles the heat flux in the throat region can be of the order 
of 20 Btu/sec/in. 2 . A heat flux of as much as 100 Btu/ 
sec /in. 2 are predicted as possible in advanced chemical 
rocket concepts. 


The nozzle calculations were made using an existing 
digital computer program^) which was revised to make 
use of these recommended equations for heat transfer 
and friction. Comparisons between predicted and meas- 
ured gas-side wall temperature, coolant pressure drop, 
and coolant temperature rise are made for several 
NERVA nuclear tests. Wall temperatures were not 
measured in the Phoebus 2 and Pewee nuclear tests, so 
that only the measured coolant pressure drop and tem- 
perature rise were compared with predictions. 

Pewee, NERVA, and Phoebus 2 vary a great deal in 
size and power but are in the same geometrical class of 
nozzles. These calculation procedures can and should 
be compared with experiments for other geometries 
when available. 

II. Basic Equations 

Coolant Side 

Heat transfer . - Heat transfer coefficients in the 
coolant passages are calculated using the straight tube 
correlation equation^ modified to include the effect of 
entrance configuration and curvature. 

The recommended straight tube correlation is 


The author of this paper has previously studied all 
available single-phase hydrogen heat-transfer data for 
flow through straight tubes and recommended a single 
correlation equation for a wide range of conditions. (1) 

The calculated heat -transfer coefficients were shown to 
be in very good agreement with the high heat flux ex- 
perimental data. ( 2 ) The straight tube equation has been 
modified to include the effects of curvature^ 2 ) and com- 
pared with experimental data for single curved tubes. ( 2 > 4 ) 
These results and recommended applications of the eq- 
uations have been reported. ( 2 ) 

In addition all available experimental friction co- 
efficients have been analyzed and a new correlation eq- 
uation recommended for single-phase turbulent flow 
through smooth straight tubes. (6) The straight tube eq- 
uation for friction can be modified to include the effects 
of curvature . ^ 

In this paper the best available heat -transfer and 
friction correlation equations are brought together and 
applied to the cooling passages of nuclear rocket nozzles. 
The most severe test of the ability of these equations to 
calculate heat transfer and friction in the coolant pas- 
sages of regeneratively cooled nuclear rocket nozzles is 
to compare the calculations with experimental results. 



where 



( 1 ) 


Eq. (1) modified to include entrance effects gives 



where F = 2. 3 for a 45° angle bend and 5 for a 90° angle 
bend. 


The Ito correction for the effects of curvature on 
friction coefficients^ has been applied to local heat- 
transfer coefficients with good results^ using Eq. (1) 
to give 
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is 0.05 for both the throat and knuckle region because 
unlike the heat -transfer coefficient which is affected by 
only the heated surface the friction coefficient is af- 
fected by the wetted surface of the coolant passage. The 
heat-transfer coefficient would then decrease in the 
knuckle region while the friction coefficient would in- 
crease. 


(4) 

for the convex side (knuckle). 

Eqs. (1) to (4) cover all the conditions encountered 
in the cooling passages of nuclear rocket nozzles. 


The Its correction is used for both heat transfer and 
friction even though it was originally proposed for only 
friction. How well the Ifo correction applies to local 
heat transfer and friction has not been studied thor- 
oughly due to the paucity of local experimental data re- 
ported for both heat transfer and friction coefficients. 


Fluid flow . - The range of conditions for which 
friction coefficients have been reported is not so great 
as those for heat-transfer coefficients. An investigation 
using all available data found that measured friction co- 
efficients have been reported which were as much as 
three times the values predicted by conventional meth- 
ods. ^ The correlation equation recommended 
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(5) 


where Re w is the Reynolds number based on diameter 
and with viscosity and density evaluated at the wall tem- 
perature. This is the Koo, Drew, and McAdams re- 
lation modified with (T w /T b )~°- 5 . For Re w a 3000 
Eq. (5) correlated all the friction coefficients within 
±10 percent. 


The friction pressure drop is calculated from 


Ap fr = 


2G 2 fx 

Dg 


( 8 ) 


The momentum pressure drop is calculated from 


1 1 

(P A) 2 (P A) x _ 

The total pressure drop is the sum of the friction and 
momentum pressure drops. 

Hot-Gas Side 


Ap mom gA 


(9) 


In calculating the heat transfer from the hot gas to 
the cooled nozzle wall, it is common practice to use the 
Nusselt equation 


Eq. (5) can be used at x/D as low as 3 if the ex- 
ponent of T w /T b is changed from 0. 5 to the c^ used 
for heat transfer. 

The Karman-Nikuradse relation can also be mod- 
ified with (T w /T b ) Cl to give 


= [4. 0 tog (RewVfj - 0 • 40] ^ ^ (6) 

Assuming that the effect of T w /T b on the friction 
coefficient is the same for a rough tube as it is for a 
smooth tube, the equation for predicting friction coeffi- 
cients in a rough tube 7 ^ becomes 



where e is the relative roughness of the tube. 

The Ito correction for the effect of curvature 
should be used with both Eqs. (6) and (7). For friction 
coefficient calculations the exponent of 


Nu f = CgRej? ’ 3 Prj? ' 4 (10) 

and vary the C g as a function of axial location or area 
ratio. The value of Cg can either be measured 1 ' or 
calculated, e.g. Bartz boundary layer calculations. 
These Cg values are usually adjusted after rocket tests 
to give calculated gas-side waU temperatures equal to 
measured values. In cooled nozzles the C„ thus be- 
comes a function of coolant side heat transfer equations. 
The theoretical values of Cg calculated by the Aerojet 
General Corp. for NERVA and Phoebus 2 nuclear rims 
are shown in Fig. 1. 

The value of Cg increases in the nozzle chamber 
from knuckle to reactor core face but the area ratio 
remains constant which is indicated by the vertical line 
in Fig. 1(a). The area ratio of the nozzle chamber is 
not the same for NERVA, Phoebus 2, and Pewee. 

Figure 1(b) shows the theoretical and adjusted Cg 
curves for the NERVA nozzle. The adjusted value of 
Cg is about 2. 5 times the theoretical value in the 
knuckle region. The need for such a high Cg in this 
region probably results from the lack of previous 
coolant side correlation equations to take into account 
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There is a factor of 3 between the two friction coeffi- 
cients. 


the decrease in heat-transfer due to the convex curva- 
ture. This effect of curvature appears in Eq. (4). The 
Cg values shown for the Pewee nozzle were experimen- 
tally determined by Los Alamos Scientific Laboratory. 

The data points were reported by Schacht, Quentmeyer, 
and Jones. ^ 

In Eq. (10) the Reynolds number is based on nozzle 
diameter. It would be more logical to treat the constant 
diameter chamber section as a flat plate and base the 
Reynolds number on axial distance along the wall were it 
not for the introduction of cold hydrogen for peripheral 
film cooling of that section . 

in. Nozzle Calculations 

The nozzle calculations were made using the digital 
computer program of Rohde, Duseha, and Derderian^ 
modified to incorporate the equations recommended in 
this paper. 

Coolant Side 

Heat transfer . - A Nusselt type correlation equation 
with the physical properties evaluated at the film tem- 
perature is part of the computer program. This eq- 
uation was replaced by Eq. (1) with the effects of en- 
trance configuration and curvature being applied in the 
proper regions as shown in Fig. 2. 

Fluid flow . - The friction coefficients are calculated 
in the program using the Karman-Nikuradse relation for 
smooth tubes and the rough tube equation which was mod- 
ified to become Eq. (7). The relative roughness e, of 
the coolant passages is as follows: NERVA: e = 5 pin. , 
Phoebus 2: e = 10 pin. , Pewee: e * 20 pin. , and Eq. (7) 
was used in the nozzle calculations as shown in Fig. 2. 

The surface temperature of the coolant passage 
varies as much as 1000° R from the crown (the semi- 
circular portion of the U-tube) to the pressure shell 
(see Fig. 3). This temperature variation makes the cal- 
culation of a single friction coefficient for the complete 
wetted perimeter very difficult. Calculations of pressure 
drop through the coolant passages are often 60 percent 
lower than the measured drop. The following results 
indicate the variation in friction coefficients which can 
result from two different assumptions for a section of 
coolant passage well downstream of the entrance region 
and upstream of the throat curvature in the NERVA 
nozzle using Eq. (7). 

Assumption 1 - The cooling passage wall tempera- 
ture is constant and equal to the hydrogen bulk temper- 
ature T w = T b . • . T w /T b = 1, Re,^ = Ret, = 2.7X10 7 
f = 0.00267. 

Assumption 2 - The cooling passage wall tempera- 
ture is constant and equal the crown temperature. T w = 
1200° R, T b = 80° R, T w /T b = 15, Re w = 8. 9X10 4 f = 
0.0009. 


The method used herein was to calculate one friction 
coefficient for the bulk temperature (f b ) and another for 
the hot crown temperature (f c ) for each station. Then an 
effective friction coefficient was defined as 



where f e ff is the effective friction coefficient, Sj, and 
S c are the surface areas at bulk and crown temperatures 
respectively as shown in Fig. 3. 

Values of (Sc/Sfc, + Sc) and (Sb/Sb + Sq) were calcu- 
lated at each calculation station along the coolant passage 
and (S c /Sb + S c ) varied from 0.40 at the coolant entrance 
to 0 . 22 at the nozzle throat and 0 . 60 in the chamber of 
the NERVA nozzle. 

Eq. (11) should be as applicable to a rectangular (or 
any other shape) coolant passage as to the U-tube pas- 
sage. 

The momentum pressure drop equation in the com- 
puter program remained unchanged. 


The heat transfer coefficients on the hot gas side 
were calculated with the Nusselt type equation in the 
original program. The values for Cg are input data 
and therefore the effects of the various Cg's can be ob- 
served. 

Hot gas static temperature and pressure calculations 
for hydrogen at equilibrium conditions for assigned 
chamber pressure, chamber temperature, and nozzle 
area ratios were made using the computer program of 
Ref. 10. 

The effect of peripheral film cooling in the chamber 
was calculated using the method of Hatch and Papell. 

In the NERVA and Phoebus 2 nozzles the film coolant 
temperature was about 1200° R and the flow rate was 
about 1 percent of the hot gas flow rate. The coolant 
was assumed to effect the coolant passage wall temper- 
ature for a distance of 6 in. from the reactor core. The 
Pewee nozzle had a coolant flow rate of about 10 percent 
of the hot gas flow rate and a temperature as low as 
200° R. For the Pewee calculations the coolant was 
assumed to have some effect all the way to the knuckle 
region. 

IV. Discussion of Results 

The only experimental wall temperatures available 
from nuclear tests are for the NERVA nozzles NRX-A3, 
EST, A5, and A6. Coolant pressure drop and tempera- 
ture rise were also measured for these tests as well as 
the peripheral coolant temperatures and flow rates. 

This information made the NRX series suitable for the 
test of the modified computer program. 


Hot-Gas Side 
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Gas-side wall temperatures predicted using the 
AGC theoretical Cg, the AGC adjusted Cg, and a con- 
stant Cg. of 0.026 are compared to experimental gas- 
side wall temperatures for NRX-A6 in Fig. 4. The ex- 
perimental wall temperatures resulted from an exam- 
ination of braze alloy patches at various locations in the 
nozzle. The open and solid symbols are at different 
angular locations (105° apart) in the nozzle. The solid 
lines represent the wall temperatures calculated with 
and without film cooling using the theoretical values of 
Cg for the NERVA nozzle geometry and is in very good 
agreement with the measured values for most of the 
length. The wall temperatures calculated with and with- 
out film cooling using the adjusted values of Cg were 
much higher than the measured wall temperatures in the 
chamber and knuckle regions. These higher calculated 
wall temperatures resulted from a combination of high 
values of Cg in the adjusted curve (see Fig. 1(b)) and a 
decrease in the coolant heat-transfer coefficient due to 
the convex curvature effect in Eq. (4). It appears that 
the effect of convex curvature has been included twice, 
once on the heat-gas side and once on the coolant side. 

It is very difficult to compare the predicted and 
measured wall temperatures near the core because the 
cold hydrogen used for peripheral film cooling is intro- 
duced at this point. Predicted wall temperatures are 
shown with and without peripheral film cooling. Ex- 
perimental studies have shown that the inclusion of a 
simulated reactor core can increase the heat-transfer 
coefficients in a nozzle chamber. 

The dashed line in Fig. 4 represents the wall tem- 
perature calculated using a constant Cg of 0.026, the 
simplified method represents the wall temperatures 
fairly well. The predicted wall temperatures near the 
core were lower than the measured temperatures even 
though film cooling (the film, coolant temperature is 
about 500° R higher than the wall temperature) was not 
included in these calculations. However, in a nozzle 
without peripheral film cooling, the use of a constant 
Cg would predict wall temperatures which could be 
much lower than the actual temperatures near the reac- 
tor core. 

For the NRX-A6 test the pressure drop calculated 
using Eq. (11) was 139 psi compared to the measured 
value of 141 psi. The calculated temperature rise is 
95° R which is a good agreement with the measured 
rise of 93° R. If the coolant passage wall temperature 
is assumed to be constant and equal to the bulk temper- 
ature, the calculated pressure drop is 153 psi, about 
8 . 5 percent higher than the measured value and still in 
good agreement with measured values. The pressure 
drop and temperature rise for the various NRX test 
nozzles are shown in Table 1. The variation of coolant 
temperature and pressure along the length of the nozzle 
is shown in Fig. 5. 

The computer program was also used to calculate 
wall temperatures and coolant conditions for the Pewee 
and Phoebus-2 nuclear tests. The coolant pressure 
drop and temperature rise for Phoebus 2 and Pewee 
tests are also shown in Table 1. 


For the Pewee nozzle the agreement between calcu- 
lated and measured temperatures rise and pressure drop 
was usually better for the full power tests than for other 
tests. 

The film coolant flow rate was about 10 percent of 
the hot gas flow rate for the Pewee nozzle compared to 
about 1 percent for NERVA and Phoebus 2 nozzles. In 
addition the coolant temperature was much lower for 
Pewee than for NERVA and Phoebus 2. These facts 
might tend to lower the confidence level for any Pewee 
wall temperature calculations to some degree. 

The fact that calculated pressure drops are in almost 
all tests about 10 percent lower than the measured value 
indicates that possibly the method of weighting the fric- 
tion coefficients on the surface area with which it is in 
contact over corrects the effective coefficients in 
Eq. (11). Another possibility is that the effect of wall to 
bulk temperature ratio which has been verified for 
smooth tubes with temperature ratios up to 7. 35 cannot 
be extrapolated to rough passages with temperature 
ratios of 15. 

V. Conclusions 

It appears from the results of this study that the use 
of the best available prediction equations for heat trans- 
fer and friction coefficients for turbulent flow through 
tubes can, with modifications for entrance and curvature 
effects, be applied with good results to the coolant 
passages of nuclear rocket nozzles. These equations 
should also be useful for any type of rocket nozzle in 
this geometry class. 

It should be concluded from the results presented in 
this paper that existing variable property heat transfer 
and friction equations can be used to predict coolant 
pressure and temperatures and nozzle wall temperatures 
with reasonable accuracy. These results form a basis 
with which future calculational methods can be compared, 
and hopefully improved. 

VI. Symbols 

A cross-sectional area 

Cg coefficient in heat -transfer equation 

c-^ exponent of T w /Tjj 

D diameter 

e relative roughness of surface 

F entrance effect coefficient 

f friction coefficient 

G mass flow rate per unit cross-sectional area 

g gravitational conversion factor 

Nu Nusselt number 
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Pr 


Prandtl number 


Ap total static pressure drop 

Apf r friction pressure drop 

Apmom momentum pressure drop 
R radius of curvature 

Re Reynolds number 

r inside radius of passage 

S surface area 

T temperature 

w mass weight of flow 

x linear distance from entrance 

p density 

Subscripts: 
b bulk 

o crown 

g gas 

w wall 
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TABLE 1 CALCULATED AND MEASURED TEMPERATURE RISE AND PRESSURE 
DROP ACROSS COOLANT PASSAGES FOR NUCLEAR TESTS 

NERVA 



EP-IV 

Anointed 


Pressure drop 

Temperature rise 






time 

Calculated 

AP C , 

psi 

Measured 

Ap m , 

psi 

Ap - Ap 

— £ -x 100, 

A Pm 

percent 

Calculated 

at c , 

OR 

Range of 
measured 

AT m , 

°R 

1 

■ 27 

29 

-6.9 

26 

13 to 36 

2 

60 

66 

-9.1 

31 

28 to 41 

3 

132 

151 

-12.6 

55 

38 to 62 

4 

182 

204 

-10.8 

52 

33 to 60 

5 

211 

255 

-17.3 

78 

57 to 88 

6 

246 

290 

-15.2 

83 

60 to 93 


Pewee 


EPIII 

Pressure drop 

Temperature rise 


Calculated 

Measured 

Ap c ' Ap m 

Calculated 

Measured 

AT C - AT m 


Ap c > 

psi 

Ap m > 

psi 

xiuu, 

Ap m 

percent 

AT C , 

OR 

AT m , 

°R 

“ X 1UU • 

AT m 

percent 

37 

59 

62 

-4.8 

62 

56 

10.7 

FP1A 

117 

126 

-7.1 

94 

91 

3.3 

B 

108 

118 

-8.5 

93 

92 

1.1 

C 

113 

124 

-8.9 

94 

93 

1.1 

D 

113 

124 

-8.9 

94 

93 

1.1 

341 

31 

37 

-16.2 

23 

24 

-4.2 

FP2A 

113 

123 

-8.1 

94 

93 

1.1 

B 

113 

125 

-9.6 

94 

96 

-2.1 

C 

113 

125 

-9.6 

94 

96 

-2.1 

D 

113 

125 

-9.6 

94 

97 

-3.1 

342 

31 

37 

-16.2 

23 

25 

-8.0 
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NERVA (THEORETICAL-AGC) 

PHOEBUS 2 (THEORETICAL-AGC) 

PEWEE (MEASURED-LASL) 

o COPPER HEAT-SINK NOZZLE (LeRC) 




Figure 1. - Variation of various correlation coefficients for 
hot -gas side heat transfer as a function of area ratio. 
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Figure 2. - Five regions of the coolant passage and equations applicable to each region. 
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Figure 3. - Typical U-tube coolant passage of a 
nuclear rocket nozzle. 
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Figure 4. - Comparison of calculated gas-side wall temperatures 
with measured values for NRX-A6 nuclear test, NERVA nozzle. 
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Figure 5. - Variation of coolant temperature and pressure 
with axial location. NRX-A6 nuclear test, NERVA nozzle. 
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